Quantum Phase Transitions and Heat Capacity in a two-atoms Bose-Hubbard Model 
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We show that a two-atoms Bose-Hubbard model exhibits three different phases in the behavior of 
thermal entanglement in its parameter space. These phases are demonstrated to be traceable back 
to the existence of quantum phase transitions in the same system. Significant similarities between 
the behaviors of thermal entanglement and heat capacity in the parameter space are brought to 
light thus allowing to interpret the occurrence and the meaning of all these three phases. 
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A Quantum Phase Transition (QPT) is a phenomenon 
[l| consisting on a qualitative change of the ground state 
of a quantum system as one or more of its physical pa- 
rameters are varied [2t51- These phase transitions are 
genuinely quantum mechanical effects directly stemming 
from quantum fluctuations. Differently from thermal 
phase transitions indeed, they occur at zero temperature 
and manifest themselves through the non-analytical be- 
havior of the ground state energy at a transition point [6|- 
Q in the parameter space of the physical system, around 
which some selected physical properties of the system 
experience significant, sometimes sharp, modifications 
[3, [lO] ■ For instance recent investigations on some sys- 
tems undergoing quantum phase transitions have demon- 
strated that the transition is related to a marked change 
in the entanglement [ll[ . Depending on the physical sys- 
tem considered, entanglement could have a peak, show 
discontinuous behavior or also manifest diverging deriva- 
tives at the critical points [i2|,[l3 ■ It is worth noting how- 
ever that, generally speaking, such transitions are accom- 
panied by correlations whose nature is not clear enough 
[lil - llq . Since the degree of entanglement getting estab- 
lished between different components of a system cannot 
be measured in laboratory, the ability of finding easily 
measurable observables that can be used as entanglement 
witnesses represents an important challenge. This paper 
investigates on a simple but not trivial system consisting 
of few atoms in an optical lattice, bringing to light the 
existence of QPTs and their connection with the ther- 
mal entanglement characterizing the system. After that 
we find our main result, that is a link between negativ- 
ity and heat capacity behaviors in the parameter space, 
leading to a deeper understanding of the origin and the 
occurrence of a QPT. 

This work is organized as follows: in Section II we discuss 
and analyze the Hamiltonian model describing our phys- 
ical system evaluating as well all quantities of interest. 
The main results of this paper are presented in Section 
HI. Comments and conclusions are drawn in Section IV. 
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describing spin-1 atoms in an optical lattice, has at- 
tracted much attention |18l - [2l| thanks to the renewed 
interest in spin gases stemming from new experimental 
techniques allowing a detailed control over multi-particle 
states [22|. In eq. ([!]) a]^ is the creation operator of 
an atom in the i — th lattice site with spin projection 
over a fixed quantization axis equal to cr = (—1,0,1), 
n.i = ^^ '^L'^io' i^ ^^^ total number of atoms on site i 
and SJq( is the total spin operator of the i—th lattice site. 
The first term of Hbh describes spin symmetric tun- 
neling between nearest neighboring sites, the parameter 
t being the tunneling amplitude and the X]<ji> being 
extended over all pairs of neighboring sites. The second 
term describes the Hubbard repulsion between atoms, 
Uq being the contact scattering amplitude, whereas the 
term proportional to U2 differentiates scattering channels 
with Stot equals to or 2, since scattering lengths for 
these two channels are not the same [2j|. Finally e is the 
atom self energy which may be neglected in what follows. 



Let us assume that two different atoms are described 
by states centered in different lattice sites and that the 
temperature is low enough to consider the probability 
of finding an atom outside the site on which its state is 
centered almost zero. 

Under this hypothesis the tunneling term can be 
treated as a perturbation and as a consequence, always at 
low temperature, the tunneling term in Hamiltonian ([T]) 
can be replaced by an effective next-neighbors interaction 
[24J . Tunneling processes induce effective pairwise inter- 
actions between atoms on neighboring sites and the tun- 
neling term in Hamiltonian ([1} , perturbatively expanded 
up to second order in parameter t, may be substituted 



with the foUowing efFcctive one 
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in which an external magnetic field w has been added 
to the system described above. Here Si represents the 
spin operator of the particle in the z-th site and J = 
J2i Sj. The effective coupling constants Kq, Ki, K^ are 
related to the parameters t, Uq and U2 appearing in the 
microscopic Hamiltonian ((U as follows 
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and satisfy the simple equation Kq ~ Ki — K2. In the 
deep Mott-insulator phase, where ([2|) holds, the tunnel- 
ing amplitude t is of the order of some kHz but can be 
as well as small as a few Hz |25l |. 

We remark that the term ^ is an effective one describ- 
ing the physics of hopping between lattice sites in terms 
of interactions between spin operators of the atoms. Two 
interaction terms appear, a linear and a quadratic one, 
and their origin can be qualitatively understood. Let us 
indeed consider atoms whose states are centered in dif- 
ferent sites. After the hopping of one of these atoms 
towards one of its neighboring sites there will be a po- 
tential well containing two particles, thus allowing them 
to interact by s-wave scattering. On the other hand two 
spin-1 particles scattering inside a potential well may re- 
sult in two different final states with total spin equal to 
or 2 (states with total spin equal to 1 are forbidden for 
symmetry reasons). Thus, in order to properly describe 
these two different situations in terms of spin-spin inter- 
action, one needs to express the projectors Pq and P2 
{Ps being the operator projecting the state of two atoms 
into a final state with total spin equal to s) in terms of 
spin operators of the particles involved. It is possible 
to show [23| that this requires quadratic interactions to- 
gether with linear ones. In this paper we concentrate on 
the simple, but not trivial, situation in which only two 
atoms arc considered. With this further approximation 
one can look at iJ^ as the full Hamiltonian of the sys- 
tem. 

Measuring the energy in units of i, H^^ in the case of 
two atoms only can be cast in the following form 
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where 



i/ = cuJ, +t(Si- 82)+ 7(81-82)' +r/ (6) 



and / is the identity operator in the 9-dimcnsional Hilbcrt 
space of the system. We do not neglect the term propor- 
tional to identity since we want to investigate the system 
in its parameters space. 

It is worth noting that investigations on few atoms mod- 
els J26l - l28l | are interesting in their own since under par- 
ticular assumptions they describe small subsystems in 
larger lattices providing hints to understand their physi- 
cal properties [23]. On the other hand considering few 
particles systems allows us to exactly diagonalize the 
Hamiltonian and to focus on the dynamically generated 
pair wise entanglement. 



A. Thermal entanglement 



To diagonalize the Hamiltonian given by eq. 
convenient to write it in the form 



it is 



H =LoJ, + §(J2 - Si - si) + |(J2 - si - Slf + 
+rl = ujJz + §(J' - 4/) + i(J' - 4/)2 + rl 

(8) 

where J = Si + S2. H is then diagonal in the coupled 
basis \jM) of common cigenstates of J^ and J^ with as- 
sociated eigenvalues equal to j{j + l) and M respectively. 
The energy eigenvalues are 

E,M=^M+^{j{j + l)-2) + l[ij{j + l)-4f-4] (9) 

Let us suppose that our system is in a thermal state at 
T temperature. The density matrix describing the two 
spin 1 atoms in the optical lattice, can be thus written 
as 



^ z 



(10) 



where /? = l/T {kg = 1) and Z = Tr{e-P"). 
Z can be exactly given in closed form as follows 



-/3t 



2 cosh /3r ( 1 -f- 2 cosh /3a;) -|- 



+2e-'5^ cosh 2/3w + e-'^(37-2r) 



(11) 



In order to quantify thermal entanglement in the state 
(IIOJ w e use the well known negativity function A^ defined 
as 1301 



^ = ^(Ei^'i-i 



(12) 



T ^ Kl/t 

7 = K2/t 

LO = U)/t 

r = T — ^ 
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where A^ are the eigenvalues of the matrix a, partial 
transpose of p with respect to one of the two spins. 
Representing p with respect to the ordered factorized 

basis || - 11), |1 - 1), I - 10), |01), |0 - 1), |10), | - 1 - 



1), 1 00), 1 11) >, and transposing with respect to the spin 
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M± = -ie-'3(^±'^)sinh(/3T) 

P± = ;i-e-'^(^±'^) cosher 

R^ ^ _^e-/3-(e-/^- ± 3e^- + 2e-^(3^-2-)) 
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Six of the nine eigenvalues can immediately be derived 
in the form 



•^3 — ^5—2 
-^4 = -^6 = 
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while the last three eigenvalues can be obtained solving 
the secular equation associated to the 3x3 block 
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Since the algebraic expressions of these last three eigen- 
values of a do not exhibit features deserving special at- 
tention we do not give them explicitly. 
The knowledge of the eigenvalues Aj of cr, as functions 
of both temperature and the three model parameters, 
allows us to investigate, at a given temperature, the be- 
havior of negativity in the parameter space. Recently 
the role of the nonlinear effective coupling between the 
two atoms on the thermal entanglement has been qual- 
itatively investigated. In detail it has been found that 
the quadratic interaction term favors the thermal entan- 
glement that on the other hand exists only for absolute 
values of the quadratic parameter 7 < larger than a 
critical value 7c [3ll | . In presence of an external magnetic 
field (w ^ 0) the system shows two different phases, one 
correspondent to A^ = 1 and the other one to iV = 0, 
and the value of 7c at which the phase transition takes 
place is dependent on w. These results have anyway been 



obtained neglecting from the very beginning the linear in- 
teraction term whose intensity is measured by t in Hamil- 
tonian model ([6|). In what follows we investigate on the 
role played by the linear term. In particular we wonder 
whether there is competition between the linear and the 
quadratic effective interaction terms or if they cooperate 
in establishing entanglement. 



III. RESULTS 

The negativity A'^ against the linear interaction param- 
eter T, for Lo and 7 fixed, is reported in Fig[T]in correspon- 
dence to three values of temperature T , namely T = 0.05 
T = 0.6 and T = 1. The addition of the term propor- 
tional to T to the model of Ref 



3l| is at the origin of an 
intermediate phase (plateau of the negativity, iV = i) for 
Lo and 7 equal to 1, when r runs around 2. N undergoes 
indeed two phase transitions (vertical lines) differently 
from what happens when r = where only one phase 
transition is present (cf. Fig.l of Ref. [3l|) It is easy 




FIG. 1: System negativity plotted against r when 7 = a; = 1. 
Energy is measured in units of t and fcs = 1 

to interpret the occurrence of the intermediate phase at 
iV = i as due to the magnetic field assisted competi- 
tion between the quadratic term and the linear one when 
both 7 and r are positive. To appreciate this point it is 
enough to note that the two-atom system, when t > 0, 
tends to minimize its total energy assuming j = (maxi- 
mally entangled state), whilst it tends to maximize both 
j and \M\ (factorized state) when 7 > in presence of a 
magnetic field. 

Figure 1 shows that the two phase transitions are rel- 
atively sharper in correspondence to the smallest value 
of the temperature T . Investigating on the behavior of 
the negativity N when the temperature T goes toward 
zero, it is possible to demonstrate that the smooth phase 
transitions shown in FigT, become indeed jumps. At the 
light of recent papers [3, Is^, we wonder whether this 



behavior reflects the existence of QPTs for the system, 
in this way taking part to the topical debate about the 
connection which might exist between entanglement and 
QPTs. It has been proved for example that in a class 
of spin systems QPTs are signalled by criticalities in the 
concurrence function adopted to measure bipartite en- 
tanglement. To this end we analyze the ground-state 
energy of the Hamiltonian model given by equation ^ 
against the parameter r. Exploiting eq. ([9|) it is possible 
to demonstrate that, in correspondence to the two tran- 
sition points of the function N, the ground state energy 
undergoes two level crossings, which are necessary and 
sufficient conditions for the existence of QPTs . We may 
thus conclude that the results we have obtained for the 
negativity N in correspondence to low, but not zero, val- 
ues of T, stem from the existence of critical points in the 
space of the parameters characterizing our system and 
thus are signatures of the presence of QPTs. 

It is in addition important to stress that the first tran- 
sition (the one from A^ = to A^ = 0.5) occurs when 
2t = Lo. However a positive value of 7 partially destroys 
the quantum correlations between the two spins not al- 
lowing the system to reach a maximally entangled state. 
Finally when r > uj+3"f the linear interaction can prevail 
on both the quadratic interaction and the external mag- 
netic field such that the system can move to a maximally 
entangled condition. This behavior does not depend on 
our choice 7 = w = 1 but manifests itself qualitatively in 
the same way whatever the two parameters arc fixed at. 

Recent papers [3J, [SJ] have demonstrated that the be- 
havior of some thermodynamical quantities, such as for 
example the internal energy, the magnetic susceptibility 
or the heat capacity, can reveal entanglement between 
the microscopic constituents of a macroscopic sample, 
suggesting an intriguing approach of experimental inter- 
est to witness entanglement. We are thus stimulated to 
search a thermodynamical quantity of transparent phys- 
ical meaning reflecting the features exhibited by the neg- 
ativity function in the system parameter space as they 
appear in Fig.l. To this end we analyze the heat capac- 
ity , which is deflned as 
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where U = T'^ 



)\nZ 
dT , 



is the mean energy of the system in 



the thermal state ([TU]). 
Figure 2 reports the dependence of Cy on t, straight- 
forwardly got exploiting the closed form of Z given by 
eq. (fTTj) . together with the negativity. The flxed value 
of temperature T in these plots is 0.6 and has been cho- 
sen taking into account the fact that a temperature close 
enough to zero would have given a flat almost zero Cy, 
while an high enough temperature would have destroyed 
thermal entanglement indeed making unsuccessful our at- 
tempt to compare negativity and heat capacity. 
In order to realize a wider intermediate phase we flx 7 
greater than one, used in FigHJ so that the transition to 
the upper phase will occur at a greater value of r. Fig. 2 




FIG. 2: Heat capacity and system negativity plotted against 
r when T — 0.6, 7 = 7 and uj = 1 



describes our main result demonstrating the existence of 
a strong similarity between the behaviors of the heat ca- 
pacity of the system and the one of negativity against 
T. As the negativity function, also the heat capacity is 
characterized by three phases each one correspondent to 
a well defined value. Moreover the higher the entangle- 
ment is, the lower is heat capacity and viceversa. It is 
worth noting that such a result is a general property that 
is independent on the region of the parameter space (r, 
7, w) here analyzed for simplicity. 



IV. DISCUSSION AND CONCLUSIONS 

Let us begin by providing reasonable arguments in fa- 
vor of our choice of performing a comparison between 
negativity and heat capacity. To this end we report 
in Fig. 3 the structure of the low-lying energy states of 
Hamiltonian ^. 

The presence of points (A and B) where clearly the 
derivative of the ground state energy undergoes a finite 
jump explains why ta and tb play the roles of quan- 
tum critical points in the behavior of negativity. It is 
worth noting that Fig. 3 shows in addition the occurrence 
of level crossings other than A and B around ta and tb ■ 
As a consequence, the gaps between the low-lying en- 
ergy levels exhibit structural changes going to zero and 
then rapidly increasing again. This fact suggests to test 
such a behavior through the heat capacity, which is by 
definition very sensitive to the structure of the low-lying 
energy levels at low temperature. The bouncing behav- 
ior of heat capacity in the vicinity of quantum critical 
points TA and tb as reported in Fig. 2 is then easily inter- 
preted in terms of these oscillations in the energy gaps 
between the accessible states of the system. The closer 
these levels are, the easier is for the system itself to ex- 
change energy with the bath and the higher heat capacity 
is. Far from the crossing points, where the ground state 
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FIG. 3: Low-lying energy levels of H (in units of t) versus r 
when 7 = 7 and oj = 1. Two crossing points, A and B in 
the plot, are clearly visible in correspondence to the values 
TA — \ and TB — 22 



energy is analytical, these gaps are constants (some level 
crossing may occur between high excited states, but these 
ones can not affect the thermodynamics of our system as 
long as the temperature is low enough) and heat capac- 
ity is thus independent on r: when t < ta there exist 
two levels, one of which plays the role of ground state, 
whose gap is constant and smaller than the mean ther- 
mal energy at T = 0.6. Heat capacity is thus constant 
and different from zero. The same can be said for the 
heat capacity when ta < t < tb, while for t > tb the 
gap between ground state and first excited level grows 
unbounded leading to zero heat capacity. 
Another interesting feature related to the comparison be- 
tween negativity and heat capacity is the observation 
that higher negativity is associated to lower thermal ca- 
pacity and viceversa. It is possible to interpret this fea- 
ture: in fact, the higher entanglement is, the stronger the 
two atoms are correlated; these strong correlations then 
force the system to stay in its ground state not allow- 
ing any transition to excited states after thermal fluc- 



tuations since the atoms can not correlate to the bath 
for monogamy reasons [35|, thus resulting in a zero heat 
capacity. On the contrary when the two spins arc not cor- 
related (factorized ground state), nothing prevents such 
a transition to happen and heat capacity reaches its max- 
imal value since now the system is able to exchange ther- 
mal energy with the bath. Here we freely speak about 
ground state properties even if our system is in a mixed 
thermal state. This is legitimated by the temperature- 
independence of the existence of a quantum critical point 
which in turn only depends on the energy spectrum prop- 
erties as seen in the parameter space. 
Moreover the temperature range we are working in (from 
T ~nK up to T ~ lOOnK) is such that quantum critical 
effects are still clearly evident, showing how almost all 
thermal and statistical features indeed directly originate 
from lowest energy levels properties. 
In conclusion, we wish to give some final remarks. 
The novel result of this paper is reported in Fig. 2 where 
we establish a direct conceptual link, in the parameter 
space, between negativity and heat capacity, giving in 
addition a transparent physical interpretation of it. We 
claim indeed that this paper for the first time successfully 
seeks signatures, at a thermodynamical level, of the ther- 
mal entanglement get established in a two-atom Bose- 
Hubbard model at low temperature, finding a surprising 
systematic correspondence at the critical points between 
quantum phase transitions and a peculiar oscillating be- 
havior of the heat capacity. 

In connection with the simplicity of the Hamiltonian 
model on which our exact results are based, we wish to 
emphasize that such kind of models has recently gained 
interest for example in the contest of the study of tun- 
neling phenomena |28|, mostly connected to Josephson 
tunneling between spin condensates [36|. It is moreover 
appropriate to underline that our physical analysis, as re- 
ported in this section, relics on arguments which might be 
useful in the interpretation of analogous physical proper- 
ties related with more complex physical scenarios (many- 
atom Bose-Hubbard models, central spin systems, spin 
chains, etc). 
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